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Spin qubits have been successfully realized in electrostatically defined, lateral few- 
electron quantum dot circuits [l]-[4]. Qubit readout typically involves spin to charge 
information conversion, followed by a charge measurement made using a nearby bi- 
ased quantum point contact [H [5l |6j. It is critical to understand the back-action 
disturbances resulting from such a measurement approach [71 |8]. Previous studies 
have indicated that quantum point contact detectors emit phonons which are then 
absorbed by nearby qubits [9l- ll3| . We report here the observation of a pronounced 
back-action effect in multiple dot circuits where the absorption of detector-generated 
phonons is strongly modified by a quantum interference eff*ect, and show that the phe- 
nomenon is well described by a theory incorporating both the quantum point contact 
and coherent phonon absorption. Our combined experimental and theoretical results 
suggest strategies to suppress back-action during the qubit readout procedure. 

The back-action process considered in this paper involves deleterious inelastic tunneling events 
between two adjacent dots in a serial double or triple quantum dot (DQD, TQD). The energy 
difference A between the initial and final electronic dot states is provided by the absorption of a 
non-equilibrium acoustic phonon, which itself is generated by the quantum point contact (QPC) 
detector [12j. Such an absorption process between adjacent dots is constrained by the energy 
conservation condition A = ^l^^ph ('^ph is the sound velocity, q the phonon wavevector). More 
subtly, it is also sensitive to the difference in phase, A(/? = d- of the associated phonon wave 
between the two dot positions, with d being the vector connecting the two dot centers [HI [T5] . 
This q (and hence A) dependent phase difference controls the matrix element for phonon-absorption 
since it determines whether the electron-phonon couplings in each of the two individual dots add 
constructively or destructively (see Fig.[T]) [16j. The result is an oscillatory probability for inelastic 
electron-transfer events involving phonon-absorption, with constructive interference occurring when 
A(/) = (2n + 1)ti (where n is an integer). 

Data showing a pronounced back-action effect are shown in Fig. [2^, which displays the stability 
diagram measured in charge detection for a few-electron DQD without a voltage drop between 
its left and right leads. The charge configuration of the quantum dot structures influences the 
conductance of a nearby QPC because of the capacitive coupling between the dots and the QPC. 
In order to serve as a charge detector it is necessary to drive a current through the detector 
QPC which, in turn, leads to the observed detector back-action. Multiple gates fabricated 85 nm 
above a high-mobility two-dimensional electron system (2DES) are used to define two dots and 
two QPCs (Fig. I2]l). The differential transconductance d/QPc/dl/L of the biased charge detector 
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FIG. 1: Interference in quantum dot— phonon interactions, a, The back-action charge fluctuations of a QPC 
(right) used to measure quantum dots generates locally non-equilibrium phonons (yellow and red squiggly lines); 
phonons emitted in the correct direction (red) can travel from the QPC to the dots and excite them from their 
ground state. In a semiclassical picture, the displacement wave associated with an excited phonon mode will have a 
maximal effect when it is exactly out-of-phase at the two dot sites (as shown in green), as it will cause an oscillation 
in the effective energy detuning between the dots. The relative phase of the wave between the dots is Acp — q - 
where q is the phonon wavevector, and d is the vector connecting the two dot centers; constructive interference 
occurs when A(/? = (2n + l)7r with n an integer. In contrast, a minimal effect is expected when the displacement 
wave is in-phase at the two dot sites (as shown in black), b, In a fully quantum description, absorption of a single 
phonon of wavevector q can occur via either the right or left dot; the amplitudes for each process add coherently 
to determine the final excitation probability. As ground and excited states can have different electronic probability 
distributions (indicated in gray), excitation leads to a measurable change in the current through the QPC charge 
detector. The right barrier is very opaque (thick vertical line) to suppress tunneling between the right dot and the 
right lead. The blue arrow indicates the excitation while the grey arrow indicates the charge transfer between dots, 
c, Schematic showing the two interfering processes for absorption of back-action generated phonons (red squiggly 
lines) by a DQD. The relative phase between the amplitude of each process is tt + A(/? (c.f. Eq. ([T])). Note that in 
a given transition, the magnitude of q is determined by the energy splitting A between ground and excited states, 
while the direction of q is largely determined by the placement of the QPC with respect to the DQD axis. Note that 
an analogous interference effect involving photon absorption is not possible in our system, as the wavelength of a 
resonant photon would far exceed the size of the nanostructure. 
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FIG. 2: Interference in back-action in a DQD. a, Charge stability diagram of a DQD at Vqpc = —1200 /xV, 
showing the differential transconductance d/gpc/dVL (in arbitrary units) numerically derived from the dc-current 
/qpc as a function of the voltages applied to gates "L" and "R". The number of electrons occupying the left and 
right dots is shown as (A^l, N-r) for ground state configurations. At the charging lines (dark; local transconductance 
minima) the overall charge A^l+A^r changes by one electron and at the charge transfer lines (white, local maxima) one 
electron moves between the dots. A triangle of back-action contains regular stripes which correspond to oscillations 
between the metastable configuration (1,1) and the ground-state configuration (1,2). b, Trace along the red line in a. 
c, Sketch of a possible back-action-induced excitation process (1,2) (2,1) (1,1) for (1,2) being the ground state 
configuration and the right tunnel barrier being closed (thick black vertical beam), d, Scanning electron micrograph 
of the sample structure, the approximate positions of the dots are marked in blue. The crystalline direction [110] is 
marked by an arrow, e, Three charge stability diagrams for Vqpc = —700 /iV, Vqpc = —900 /xV, and Vqpc = —1400 
/iV. The blue lines enclose the total striped area which increases proportionally to Vqpc- The dashed lines in a, b, 
and e correspond to A = |eVQpc|- f? Detuning Amax corresponding to the size of the striped triangle as a function 
of Vqpc for two different series of measurements for differently tuned DQD systems (data in e belong to the blue 
circles). The error bars contain both the uncertainty of determining the voltage-to-energy conversion (see Methods 
section) as well as the error in triangle size from the measured data. The straight line denotes Amax = |eVQPc|. 
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QPC (Vqpc = — 1.2mV) is plotted as a function of control gates Vl^^r- It shows local extrema 
at the boundaries between regions of different electronic ground states, yielding dark "charging" 
and white "charge transfer" lines. Specific ground state configurations are labeled (A^l, ^r), where 
the integer Na denotes the number of electrons in dot a = L (left) and R (right). As our DQD is 
cooled to T ~ 30 mK the unmeasured DQD is expected to be in its ground state. 

Detector back-action manifests itself within a distinct triangular-shaped region of deviations 
from the ground state configuration (1,2), where a pronounced pattern of repeated, parallel stripes 
is present. It indicates an oscillating probability to find the DQD in the excited configuration (1,1). 
The excitation process sketched in Fig. ^ includes an inelastic tunneling transition (1,2) (2, 1) 
mediated by the absorption of a phonon, followed by an elastic (and therefore quick) tunneling 
process (2,1)^ (1,1). In our measurements, the tunnel barrier between the right dot and right 
lead is tuned to be almost closed (see Figs. [2]3,d). The direct transition (1, 1) (1, 2) back into the 
ground state via an elastic tunneling process from the right lead is consequently very slow and the 
excited configuration (1,1) is metastable. The associated three-level dynamics can result in average 
non-thermal occupations [T3]. In this way a metastable excited state is essential to directly observe 
detector back-action in a low bandwidth stability diagram measurement. It requires asymmetric 
dot-lead tunnel couplings in case of a DQD (cf. Supplementary Information). 

The stripe pattern constitutes the key signature of the coherent phonon-mediated back-action 
effect. It indicates that the probability to be in the excited configuration (1,1) oscillates as a 
function of the energy detuning A between the intermediate state (2,1) and the ground-state 
configuration (1,2) (see Fig. each stripe is thus parallel to the white charge transfer line 

where these states are degenerate (i.e. A = 0; marked). The striped region is bounded by a line 
A = Amax, indicating that there is a maximum energy available to excite the DQD. By seeing how 
this boundary changes with increasing Vqpc (Figs. [2^-f), we find that A^ax — |eVQPc|, consistent 
with the QPC indeed being the energy source for the initial DQD excitation. 

The geometry of the back-action regions as well as the influences of temperature and the orbital 
excitation spectrum are discussed in the Supplementary Information. In short, the remaining 
boundaries of the triangle-shaped regions of back-action correspond to energy thresholds for lead 
tunneling. The width of each stripe is largely independent of temperature; this is indicative of 
an excitation process involving electron transfers between dots, without any involvement of lead 
electrons (see Fig. ^ and Figs. [3^,f). The regular spacing of the stripe features in both DQD 
and TQD (discussed below) experiments over so many stripes eliminates the possibility that they 
are due to resonances with orbital excitations of a dot, as there is no reason to expect such a 
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FIG. 3: Interference in back-action in a TQD. a— d, Charge stability diagrams of a serial triple quantum dot 
(differential transconductance d/gpc/dKo as in Fig. [2]). The linear response of /qpc is measured while the voltage 
applied to the gate marked with is modulated by 6Vr^ = 1 mV rms at the frequency of / = 13 Hz; inset of Fig. 
[3]d: scanning electron micrograph of the sample structure. Blue dots mark approximate positions of the quantum 
dots. The number of electrons occupying the left, central, and right dots is shown as (A^l, A^c, A^r) for ground state 
configurations. The QPC is biased by Vqpc = —100, —300, —500, and — 700/iV for a-d. Additional features are 
deviations from the ground state configuration caused by detector back-action, e, Sketch of a possible back-action- 
induced process producing the striped pattern marked by the upper arrow in the (0,1,1) region in d. A phonon is 
absorbed and transfers an electron from the central to the left dot. Subsequently, the electron tunnels to the left 
lead, lowering the overall energy (of the system including dots and leads). The resulting configuration (0,0, 1) is a 
metastable excited state, f, Sketch of the back-action process in the (1,0,0) ground state region (lower arrow in d), 
which is also based on the absorption of a phonon, a short-lived intermediate, and a metastable excited state. The 
transitions in b and f are analogous to the process described in Fig. with the closed tunnel barrier replaced by 
the rightmost dot in Coulomb blockade. 

uniform level spacing; further, the energy spacing between the stripes is much smaller than would 
be expected for the average level spacing of the small dots studied here. 

To quantify the interpretation of the stripe pattern in Fig. [2] in terms of interference and QPC 
back-action, we have developed a theoretical model which describes the generation of phonons 
by the non-equilibrium QPC charge fluctuations [8J, and their coherent absorption by the DQD. 
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These fluctuations represent the fundamental back-action of the measurement- their magnitude is 
bounded from below by the rate at which information is obtained from the QPC via a Heisenberg- 
like inequality [8]. Given this, the back-action charge noise mechanism we describe must necessarily 
make a contribution to the observed oscillations. This mechanism is also consistent with the high 
visibility of the oscillations, as such visibility requires a highly localized source of hot phonons. 
While we cannot completely rule out that other, less direct back-action mechanisms contribute 
additionally (e.g. generation of hot phonons in the QPC leads), it is not clear that such mechanisms 
would also yield such high-visibility oscillations. We describe bulk acoustic phonon modes of GaAs 
interacting with both electrons in the DQD, as well as with the fluctuating charge density of the 
biased QPC via a screened piezoelectric interaction. Using Keldysh perturbation theory, we can 
calculate the DQD state in the presence of back-action (see Fig. [4] and Supplementary Information). 
The relevant part of the dots-phonon interaction (i.e. terms that can cause transitions in the dot) 
take the form: 

^i"t-|E\>(^^''"^-^^''"") (v.-+<-.-) {\9)(^\ + h.c) (1) 

Here \g) (|e)) denotes the DQD ground (excited) state, tc is the pertinent interdot tunnel matrix 
element, destroys a phonon of wavector qin branch /x, and A^*^^ is the effective matrix element 
(screened) for the interaction of phonons with a single dot. The first bracketed factor in the sum of 
Eq. ([T]) denotes the key interference of relevance: the two terms correspond to phonons interacting 
with electrons in either the left or right dot (which are centered at and respectively), see 
Fig.[lJ). 

Despite the explicit interference evident in Eq. Q, geometric averaging can still strongly sup- 
press interference oscillations in observable quantities. Simply put, while the DQD ground-excited 
energy splitting fixes the magnitude of a phonon participating in an inelastic tunneling event, it 
does not specify its direction; hence, the relative phase in the first term of Eq. ([T]) is not completely 
determined by A. This is typically the case in situations probing the emission of acoustic phonons 
by biased DQDs [16j, where interference oscillations are observed, albeit with much smaller vis- 
ibilities than seen here [TTl [T8J. In contrast, the simple geometric filtering depicted in Fig. [T^ 
suggests that this averaging need not play a role in phonon absorption, as only phonons traveling 
from the QPC to the dots contribute. This is supported by our theoretical calculations, which also 
exhibit strong oscillations for realistic parameter values, and show a pronounced enhancement of 
interference oscillations when the DQD and QPC are all collinear (see Fig. [4^-c). 

The theory is also able to capture other aspects of the experimental data: in particular, the 
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FIG. 4: Theoretical results, a, DQD stability diagram calculated using the theory described in the main text, 
which describes the emission of acoustic phonons by the non-equilibrium QPC charge fluctuations, and their coherent 
absorption by the double dot. Plotted is the derivative of the weighted average DQD charge (ur) -\-e{nL) (e — 0.4, see 
Methods section) with respect to the gate voltage Vl in the DQD experiment, as a function of the gate voltages Vl 
and Vr (see Fig. [2]); this is proportional to the measured transconductance of the QPC. Parameters are taken to be 
the same as those in the DQD experiment; in particular, the dots and QPC are taken to be collinear and aligned with 
the [110] crystallographic axis. For this orientation, one finds that the so-called "fast transverse" acoustic phonon 
mode "24 makes the dominant contribution. Pronounced stripe patterns are seen, similar to the experiment, b, A 
cut of the calculated stability diagram in a), along the indicated line. The suppression of oscillations at low values 
of ground state-excited state energy detuning A is the result of screening, and is consistent with experiment. The 
oscillations are cut off at A = Amax leVqpcl, also in agreement with experiment, c, Calculated stability diagram 
for identical parameters as a), except that the QPC is not collinear with the two dots (the QPC-dot axis is rotated 
70° from the dot-dot axis d). The result is a suppression of the back-action-induced stripe pattern, d, Calculated 
stability diagram for identical parameters as a), except that the orientation of the DQD-QPC axis is now rotated 
20° from the [110] direction. The resulting suppression of interference is due to the anisotropy of electron-phonon 
interactions: the effective electron-phonon interactions are weaker in this direction. The change in orientation also 
causes the "slow transverse" acoustic phonon mode to contribute, yielding a second, high-frequency oscillation. Insets 
in c and d show cuts through the stability diagram of the transconductance along the same line indicated in red in 
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size of the back-action triangle grows with |Vqpc|, and the lowest-energy stripes (i.e. smallest 
values of A corresponding to long phonon wavelength) are suppressed due to screening effects (see 
Fig. 4b). Using the fact that in the experiment the QPC and DQD are approximately collinear, 
the measured spacing of the interference parameter 6 A = 45/xeV in the DQD data of Fig. [2]l yields 
a DQD separation d = hvg/SA 2^ 250 nm; this is in good agreement with the separation estimated 
from SEM images (Fig. [2]l). The theory also shows that due to the anisotropy of the electron- 
phonon matrix elements A^*^^, the overall magnitude of the phonon-induced back-action is sensitive 
to the orientation of the dots-QPC axis with respect to crystallographic axes. This dependence on 
orientation is demonstrated in Fig. [4]l. More details on the theoretical treatment is provided in 
the Methods section and Supplementary Information. 

While we have focused so far on back-action in DQDs, the mechanism we describe is extremely 
general, and is in fact even more ubiquitous in systems with more than two dots. As discussed, a 
key requirement to see the effect is the existence of a long-lived metastable excited state. Such a 
situation occurs rather naturally in serial TQD structures [I9l423j . as the center dot is effectively 
decoupled from one of the leads whenever either one of the other two dots (left, right) is in Coulomb 
blockade. This directly yields a metastable excited state in which the charge of the middle dot is 
unable to relax. As a consequence deviations from the ground state configuration are often observed 
along the charging line of the center dot [20] and back-action effects occur naturally in the stability 
diagram. We study detector back-action in a TQD in Fig. [s] by successively increasing |eVQPc|. 
Already at relatively small bias |Vqpc| < 300 /xV (Figs. and b) a triangular-shaped region of 
telegraph noise is observed along the central charging line [20j. It indicates slowly fluctuating 
deviations from the ground state configuration, which can be caused by external noise or detector 
back-action [lOj. The underlying excitation processes, sketched in Figs.[3^-f, are similar as the one 
discussed above for the DQD. Indeed, the population of the right dot does not fluctuate; it takes 
the same role as the closed barrier in case of the DQD, namely to block charge exchange between 
the center dot and the right lead. Further increasing |Vqpc| to 500 /xV in (Fig. ^) reveals the 
familiar pattern of equally spaced stripes both within the (1, 0, 0) and (0, 1, 1) regions. As the bias 
is increased even more to |Vqpc| = 700 (Fig. [3]l) the striped regions expand further, revealing 
the Vqpc dependence also observed in case of the DQD (see Fig. [2]F). 

By considering experimental data on both DQD and TQD systems, we have demonstrated that 
interference can strongly affect the phonon-mediated back-action generated by a QPC in quantum 
dot circuits. Further, we have shown that this effect is well described by a basic theoretical model 
incorporating the generation of phonons by the QPC detector and their coherent absorption by 
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the dots. Our study suggests the possibihty of mitigating back-action effects by making use of 
this interference. One could, for example, endeavour to first tune the DQD/TQD to an operating 
point where destructive interference suppresses phonon absorption, and only then energize the 
QPC to make a measurement. More complex schemes which also incorporate the anisotropy of 
the electron-phonon interaction with respect to crystallographic axes could potentially yield even 
greater back-action reduction. Since the piezoelectric coupling to in-plane phonons is maximized 
in the (110) directions [24j, by aligning the QPC-DQD axis away from these directions, one could 
appreciably decrease the phonon-mediated back-action excitation discussed here (see e.g. Fig. 4a 
vs 4d). 

METHODS 

Experiment. The samples have been fabricated from GaAs / AlGaAs heterostructures con- 
taining 2DESs 100 nm (TQD) and 85 nm (DQD) below the surface, respectively. The 2DESs are 
characterized at cryogenic temperatures by the carrier densities of n^- = 2.1 x 10^^ cm~^ and 
ns = 1.9 X 10^^ cm-2 with the mobilities of = 1.72 x 10^ cm^Vs and fi = 1.19 x 10^ cmVVs for 
the TQD and DQD, respectively. Metallic gate electrodes have been fabricated on the sample sur- 
face by electron-beam lithography and standard evaporation/liftoff techniques. Negative voltages 
applied to these electrodes are used to locally deplete the 2DESs and thereby define the quantum 
dot and QPC structures. All measurements have been performed in dilution refrigerators at cryo- 
genic temperatures below 100 mK. To detect the charge configuration of the TQD the voltage of 
one gate of the TQD was slightly modulated and the detector differential transconductance was 
measured in linear response (ac set-up). In addition a constant voltage was applied across the 
QPC to enhance detector back-action. In case of the DQD only a constant voltage was apphed 
across the QPC and the direct current /qpc flowing through the QPC was measured to detect the 
charge configuration of the DQD (dc set-up). The differential transconductance d/gpc/dVL was 
then computed numerically. Both methods result in the differential transconductance and their 
interpretation is identical. 

To interpret the observed back-action in terms of the energy detuning A between different charge 
configurations an accurate conversion from gate voltages to units of energy is necessary. Such a 
linear transformation has been performed following the methods described in [25j. The conversion 
relation reads A = (^^^l ~ ^gL) + (^^^r ~ ^gR) ^' with the following set of conversion factors 
determined for the red symbols in Fig. c^^l= (54 ± 5) meV/V, a^^^ (105 ± 4) meV/V, a^^^ 
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(62 ± 4) meV/V, <^gL= (90 ± 5) meV/V. The conversion factors related to the blue symbols in 
Fig. |2|F read afj^= (65 zb 6) meV/V, a^^= (109 ± 7) meV/V, (61 ± 6) meV/V, (97 

± 7) meV/V . 

Theory. The fluctuating non-equilibrium QPC charge density operator p(f) is modelled as 
p{r) = f{r)Q, where the total charge operator Q is described by scattering theory (cf. Refs. [8,..26j). 
Note that as we are interested in a single-channel QPC, the spatial profile /(r) of the fluctuating 
QPC charge density is fixed; for simplicity, we take it to be a Gaussian of width rgpc- This 
fluctuating QPC charge density is coupled to acoustic phonons via the standard piezoelectric 
interaction (using parameters appropriate for GaAs [24j). We calculate the Keldysh Green functions 
of the acoustic phonons in the presence of this coupling to the QPC, working to leading order in 
the electron-phonon interaction, and using scattering theory to calculate the QPC Keldysh Green 
functions. We then use these "dressed" phonon Green functions to calculate the Fermi Golden rule 
excitation rate of the DQD via the coupling described in Eq. ([T]). This excitation rate is finally 
incorporated into a master equation describing the occupation probability of the three relevant 
DQD states (see Fig. In addition to the excitation rate (top panel of Fig. there is a 

rate Ffast describing the tunneling from the excited state to the metastable auxiliary state (middle 
panel of Fig. and a rate Fgiow describing the slow decay back to the true ground state (bottom 
panel of Fig. We take Ffast = 1 GHz and Fgiow = 10 kHz; in this regime of Ffast ^ ^s\ow^ 

the non-ground state population of the DQD is independent of Ffast, whereas Fsiow determines 
the overall magnitude of the interference oscillations. By using the master equation to calculate 
the stability diagram as a function of gate voltages, one can obtain the DQD charge susceptibility 
d{{T^R) + ^{TiL)) / dVL , which is proportional to the measured transconductance. The parameter 
e ^ 0.4 characterizes the QPC's asymmetric response to charge in the R versus L dot. Further 
details about the explicit form of A^*^^ (including the role of screening and dimensionality) are 
provided in the Supplementary Information. 
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SUPPLEMENTARY INFORMATION 



THE ORIGIN OF DETECTOR BACK-ACTION-INDUCED NONEQUILIBRIUM OCCU- 
PATIONS 

The stability diagrams in Figs. 2 and 3 in the main paper are modified by the absorption of 
phonons emitted from the quantum point contact (QPC) which is used as a charge detector. In this 
section we detail the microscopic tunneling processes involved in examples of detector back-action 
and discuss the general conditions under which the observed situations of nonequilibrium charge 
configuration occupation occur. 

The microscopic back-action processes involved in the triple quantum dot (TQD) are sketched 
in Figs. 3e,f. In the lower right region of back-action the ground state configuration is (1, 0, 0). As 
shown in Fig. 3f the electron in the left dot is transferred to the center dot (1, 0, 0) -> (0, 1, 0) via 
the resonant absorption of a phonon. Then, an additional electron tunnels from the left lead into 
the left dot (0, 1, 0) (1,1,0). This latter process is a very fast resonant tunneling event since the 
lead provides a continuum of occupied states below the Fermi energy. However, the configuration 
(1,1,0) is metastable since the relaxation transition (1,1,0) (1,0,0) is hindered by Coulomb 
blockade. 

In the upper left region of back-action the ground state configuration is (0,1,1). As shown in 
Fig. 3e, the electron in the center dot is transferred to the left dot (0,1,1) (1,0,1) via the 
resonant absorption of a phonon. Following this the same electron tunnels from the left dot into 
the left lead (1, 0, 1) (0, 0, 1). This is again a very fast resonant tunneling process since the lead 
provides a continuum of unoccupied states above the Fermi energy. However, the configuration 
(0,0,1) is metastable since the transition (0,0, 1) -> (0, 1, 1) is hindered by Coulomb blockade. 

Note that in both situations described above, the charge of the right dot remains fixed. Charge 
exchange between the center dot and the right lead is always hindered by Coulomb blockade. In case 
of a double quantum dot (DQD) we consequently just replace the right dot by a large tunnel barrier, 
in order to suppress charge exchange between the right dot and leads. Indeed, the microscopic back- 
action process observed in Fig. 2 is identical (except for uninvolved electrons) to the second one 
described above and can be summarized by the two-stage transition (1, 2) (2, 1) (1, 1). 

In our measurements we observe a strong nonequilibrium occupation and can even reach full 
population inversion, where the ground state remains unoccupied. The principles of thermodynam- 
ics permit such a situation, but only under two specific conditions: First, a nonequilibrium energy 
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source is required (in thermal equilibrium the ground state has the largest occupation probability 
and the occupation of all states is predetermined by the Boltzmann factor). In our case the QPC 
charge detector acts as the nonequilibrium energy source. Secondly, transitions between at least 
three states must be involved (in nonequilibrium, population inversion is impossible in a two-level 
system because the emission rate is always larger than the absorption rate due to spontaneous 
emission). This second condition emphasizes the importance of the (short lived) intermediate state 
which decays rapidly into the metastable excited state. It is the Coulomb blockade effect which 
makes this excited state metastable and enables one to observe the detector back-action. If not all 
of these requirements are fulfilled the back-action effect might still exist but would not be directly 
observed within the stability diagram of a quantum dot circuit. 

GEOMETRY OF BACK-ACTION REGIONS IN STABILITY DIAGRAMS 

In the following we discuss the geometrical shape of the regions which contain back-action and 
interference patterns (parallel stripes). Fig. 2 of the main paper shows one back-action region for 
a DQD and Fig. 3 contains two distinct regions of back-action in the stability diagram of a TQD. 
As an example we replot in Fig. [5] the charge stability diagram of the DQD for Vqpc = — 0.9mV 
already shown in Fig. 2e. It shows the differential transconductance dlqpc/dVL measured with 
a biased QPC which is capacitively coupled to the DQD (compare Fig. 2d); the black hues are 
charging lines while the white line is a charge reconfiguration line. As this measurement has 
been performed at very low temperature (T2DES ^ 35 mK), in thermal equilibrium the charge 
configuration of the DQD is well defined in the areas between these lines and corresponds to the 
ground state. The deviations from the ground state configuration (1,2) in a triangular-shaped 
region in Fig. [5] are the result of detector back-action. The microscopic processes involved in this 
nonequilibrium phenomenon and the resulting interference pattern (parallel stripes) are discussed 
in detail in Sec. [] as well as the main paper. 

Here, we specifically address the position and boundaries of the triangle. For this purpose it is 
sufficient to consider the transition (1,2) (2,1) for which an energy detuning. A, needs to be 
overcome. The tip of the triangle [on the triple point (1, 1) ^ (2, 1) ^ (1, 2)] is characterized by 
A = as in the lower right sketch in Fig.[5j The triangles base line, which lies parallel to the charge 
reconfiguration line (1, 2) ^ (2, 1), is determined by A 2^ |eVQPc| (compare upper left sketch and 
red arrow in Fig. [s]). This range in A spans the full spectrum of energy quanta that the QPC can 
emit. On its left the triangle is bounded by the charging line of the right dot, (1,1) ^ (1^2). Below 
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FIG. 5: Geometry of back-action regions. Charge stability diagram of the serial DQD (excerpt from 
Fig. 2e) measured at T2DES — 35 mK using charge detection with a QPC biased with Vqpc = — 0.9mV 
(compare Fig. 2d). Each dot is charged by one or two electrons. The triangular-shaped region containing 
back-action-induced average deviations from the ground state occupation is framed by a blue line. The 
energy detuning A between configurations (1,2) and (2, 1) along the red arrow spans the regime < A < 
k^QPcl- The four sketches, illustrating the ground state configuration at selected points just outside of 
this triangle, are used to explain the geometry of the back-action region (see text). The phonon- mediated 
transition relevant for the observed interference pattern is also indicated in one sketch (cf. main paper). Blue 
areas depict occupied states in the degenerate 2DES of the leads with Fermi energy while the vertical 
lines represent the tunnel barriers. Black horizontal lines mark the discrete chemical potentials /iL and /iR of 
the dots, the energy needed for adding the third electron, based on configuration (1, 1). The tunnel barrier 
between the right dot and the right lead (fat vertical line) is almost closed. 

this line, (1,1) is the ground state configuration and the probability for (1,2) to be occupied is low 
(see lower left sketch in Fig. §. As a consequence the transition (1, 2) (2, 1) is suppressed and 
no back-action is observed. The right boundary of the triangle is a continuation of the charging 
line of the left dot (1, 1) ^ (2, 1). Above this line the configuration (2,1) has a lower energy than 
(1,1) (compare upper right sketch in Fig. [s]). Here, the transition (1,1) (2,1) occurs rapidly 
via resonant tunneling of an electron from the left lead into the left dot; the configuration (1,1) is 
no longer metastable, so the ground state configuration (1,2) is occupied most of the time and no 
back-action is observed. 

The boundaries of the detector back-action-induced triangles in the TQD stability diagram in 
Fig. 3 can be explained similarly. 
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Influence of the electronic excitation spectrum 

The back-action-induced nonequihbrium occupation within the triangles oscillates as a function 
of energy detuning A. In the differential transconductance signal this leads to a regular pattern of 
stripes parallel to the relevant charge reconfiguration line (see, for instance, Figs. 2, 3 and 4 in the 
main paper). We interpret these occupation oscillations as an interference pattern of two competing 
phonon-absorption processes which enable an inelastic interdot electron tunneling transition. The 
period 6 A of the interference pattern corresponds to the energy of a phonon with a wavelength 
matching the distance between the two dots (see main paper). The regular spacing SA ^ SO/xeV 
of the interference stripes excludes any interpretation in terms of the electronic excitation spectra 
of the dots, which are less regular with energy spacings much larger than 50 fieY. Note that small 
variations of the period 6 A are the result of a beating of different contributing phonon modes, 
details of which will be the focus of a future paper. Let us now discuss the role of the electronic 
excitation spectrum of the dots. 

Electronic excitations in the individual dots can influence the occupation probability of nonequi- 
hbrium configurations; an example is given in Ref. [13j. In our experiments discussed here, the main 
effect of excited dot states is their influence on the refilling rate, which returns the system from the 
nonequihbrium configuration back to the ground state configuration. This influence can be directly 
seen in Fig. [6^ which plots a typical charge stability diagram in differential transconductance for 
the case when the barrier between the right dot and the lead is almost closed (compare sketches 
in Fig. ^jp). Two dark lines, parallel to the charging line of the right dot, are visible in the back- 
action-induced triangle in Fig. |6^. They mark resonances of excited states of the right dot with 
the Fermi energy. The characteristic level spacing of the electronic spectrum is roughly 120 /xeV, 
considerably larger than the phonon-induced oscillation period of 6 A ^ 50 /xeV. If the chemical 
potential of such an excited state falls below the Fermi energy (above the relevant resonance in 
Fig. [6^) an additional channel for the transition (1, 1) (1,2) via elastic co-tunneling from the 
left lead (or direct tunneling through the almost closed right barrier) is opened up. The result 
is an increase in the occupation of the ground state configuration (1,2) (leading to the dark lines 
in differential transconductance). The sketches in Fig. ^jp depict the relevant elastic channels of 
the transition (1,1) (1,2). The deviation from the triangular shape of the back-action-induced 
region at large detuning in the lowest panel of Fig. 2e in the main paper is accordingly explained 
by an excited state in the right dot which strongly enhances the transition (1,1) (1,2). 

The interference pattern in Fig. [6^ stems from transitions between the electronic ground states 
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refilling 



FIG. 6: Influence of the electronic excitation spectrum of the DQD. a, Charge stabiHty diagram 
measured with the DQD sample at Vqpc = — l.OmeV, similar to Fig. 2a. Two electronic excited states in 
the right dot cause additional dark lines of negative differential transconductance parallel to the left border 
of the striped triangle, b, Same plot as in a, with the additional lines marked in red. The upper sketch (as in 
Fig. |5| shows the phonon-mediated transition (1, 2) (2, 1), which primarily involves the electronic ground 
states of the two dots and which is immediately followed by the quick transition (2, 1) (1, 1) (compare 
Sec. [|. The other two sketches depict possible processes of the refilling transition (1, 1) (2, 1) if only the 
ground state (one channel) or in addition an excited state contributes (two channels). 

of the two dots mediated by the absorption of a phonon. Intradot transitions involving excited 
electronic states would cause characteristic phase shifts in the interference pattern, which we do 
not observe in our experiments. The two black lines in Fig.[6]cross the interference pattern without 
disturbing it. We therefore conclude that the decay of electronic excited states is fast compared to 
the interdot transition rate. 

TEMPERATURE DEPENDENCE 

The temperature dependence of a slice through a stability diagram of the TQD (compare Fig. 
3 of the main paper) is shown in Fig. [7| It includes the first charging line of the right dot (lower 
horizontal and dark line in Fig. [7| and a small number of clear interference stripes of the triangle 
in the lower right corner in Fig. 3 of the main paper. The washed out region above the clear 
interference stripes contains the first charging line of the center dot which is, however, smeared 
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FIG. 7: Temperature dependence of the interference pattern: Shown is a vertical slice through a 
stability diagram of the TQD as the ones in Fig. 3 of the main paper (but measured in a different cooldown) 
as a function of temperature for Vqpc = — 900 /iV and at constant = —1.272 V. On the left-hand side, 
the respective ground states are sketched. The dark charging line at the resonance (0,0,0) ^ (1,0,0) 
shows clear temperature induced broadening for T > 300 mK. The interference stripes above show no such 
broadening. They belong to the back-action triangle in the lower right corner in Fig. 3 of the main paper. 
Above a few clear interference stripes is the first charging line of the center dot obscured by back-action. 

out by back-action. The width of the dark charging line increases with temperature as can be 
clearly seen in Fig. [7| for T > 300 mK. This effect is caused by electrons tunneling between the left 
dot and its lead and the temperature broadened Fermi-Dirac occupation of the electronic states in 
the lead. Considerably less temperature broadening is observed for the interference stripes. Such 
a weak temperature dependence is only possible if the interference pattern is based on interdot 
transitions which are not influenced by the temperature broadening of the Fermi edge in the 
leads. This observation clearly supports our interpretation of the observed back-action in terms of 
nonequilibrium phonons driving interdot transitions. 
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THEORETICAL MODELLING 



Electron-phonon interaction 

As discussed in the main text, the back-action process of interest involves two charge config- 
urations where an extra electron is in one of either two adjacent quantum dots (e.g. in the DQD 
system described in Fig. 3, these are the charge states (2, 1) and (1, 2)). For simplicity, we present 
theory for the DQD system, and focus on the state of the extra electron. The two relevant charge 
states are thus \L) (extra electron in left dot) or \R) (extra electron in right dot). The Hamiltonian 
of the DQD takes the form 

l^dqd = e {\L){L\ - \R){R\) + {\L){R\ + \R){L\) , (2) 

where e denotes the difference in electrostatic energies between the states |L), and tc denotes 
the tunnel coupling between left and right dots. 

Restricting attention to the subspace spanned by the states |L), the Fourier transform of 
the DQD electron charge density operator at wavevector q takes the form: 

PmM = a^\L){L\^p^\R){R\^-f^\L){R\ + il^\R){L\ (3) 

where = (L|e^^'^|L), = {R\e^^'^\R) and 7^ = {L\e^^'^\R). We take the ground state wave- 
functions of the left and right dots to be Gaussians of width rg, centered at and r^, respectively. 
The system geometry is shown in Fig. [8j Similar to previous theoretical treatments of phonon 
emission by a DQD [UlETj, we neglect the coupling between the off-diagonal elements of the DQD 
density operator in Eq. ^ and phonons, as such terms only weakly modify the main effect due to 
the diagonal terms in p^. 

QD L QD R 




FIG. 8: Geometry of the double quantum dot with respect to the nearby QPC. 
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Using the expression for the charge density operator, we can now write the interaction between 
DQD electrons and phonons in the standard manner. Similar to Ref. [27j, we focus on the interac- 
tion with acoustic phonons via the piezoelectric interaction. Unlike Ref. [27j, we keep all details of 
the acoustic phonon spectrum of GaAs (i.e. anisotropic sound velocities, polarizations), as obtained 
from a standard elasticity-theory calculation [24j, and also include the effects of screening. It is 
useful to write the electron-phonon interaction in terms of the eigenstates of i^dqd which we denote 
1^) (ground state) and |e) (excited state). Keeping only those terms which can generate transitions 
between the DQD eigenstates, we obtain the interaction Hamiltonian i^int given in Eq. 1 of the 
main text. The matrix element A^^^ appearing in this equation is the effective screened matrix 
element for the piezoelectric interaction of phonons in mode with the electron density of a single 
quantum dot; the mode index fi G {L, ST, FT} refers to the longitudinal, slow transverse and fast 
transverse modes, respectively. The matrix elements take the general form \24\ I28|. 

where pM is the mass density of GaAs and V is the appropriate crystal volume element. The 
phonon frequencies are cj^^^ = Cq^^ where the corresponding sound velocities c^^^ depend on the 
wavevector direction and are calculated within elasticity theory [2l]. As discussed in the Methods 
section, the Gaussian factor e~^^^o/2 is a cutoff coming from the finite size of each dot; we take 
ro small enough (ro — 2 nm) so that this cutoff plays no role. We discuss each of the remaining 
factors in what follows. 

A4^^ is the bare (unscreened) piezoelectric coupling matrix element, and takes the standard 
form: 

-M^i = i^ehu? (qi h + h qs + qs qi et^f (5) 

where hu = 1.44V/nm is the piezoelectric constant for GaAs [29]. Here, Qi = qi/\q\ is the 
normalized component of the phonon wavevector along the crystallographic axis i G {x^y^z}^ 
while e^[q\ is the projection of the given phonon mode's polarization vector onto the z-axis. 

The form factor J^g^,^ in Eq. Q accounts for the suppression of the interaction of electrons 
with phonons having a large wavevector component normal to the plane of the 2DES. Assuming a 
standard triangular form for the transverse confining potential of the 2DES, one obtains: 



j dz\po{z)\\-"^^'\ p^{z)^e{z)^^Qe-i-a (6) 



where 9{z) is the unit step-function, a — 3.5 nm is the 2DES thickness and q± is the component 
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of the wavevector perpendicular to the 2DES plane. Here po{z) is the transverse wavefunction of 
a 2DES electron. 

Finally, the factor describes the effect of screening in the plane of the 2DES. Using a 

standard RPA approach which accounts for the two-dimensional nature of 2DES electrons [30] . 
one obtains: 

5,|„A= (^^^)', H = j dz' \po{z')\^ j dz\po{z)\\-'i^\\^-^'\ (7) 

where the effective screening radius = 5 nm is equal to half the Bohr radius of GaAs, and q\\ is 
the in-plane component of the phonon wavevector. As expected, the screening factor suppresses 
the contribution from long- wavelength phonons, and thus suppresses the back-action stripe pattern 
at small values of the energy detuning A. 

Quantum point contact charge noise spectrum 

As discussed in the main text and Methods section, the fluctuating electronic charge density 
associated with the QPC locally generates non-equilibrium acoustic phonons. The strength of 
these charge fluctuations are described by the quantum noise spectrum Sqq[uj\ of the QPC charge 
operator Q\ 

/oo 
dte-*(Q(t)Q(0)) (8) 
-OO 

This charge-fluctuation spectrum can be calculated using a standard scattering-theory approach 
to mesoscopic transport [8l|26l|31]. The relevant, negative frequency part of the spectrum (which 
describes the emission of energy by the QPC) takes the following form at low temperature : 

Here, T is the transmission of the QPC, AT (AU) is the change in QPC transmission (potential) 
resulting from changing the DQD charge state from \L) to and V^pc is the QPC bias voltage. 
Note that the magnitude of these charge fluctuations is set by the sensitivity of the QPC to the 
DQD charge state- this is a direct consequence of these charge fluctuations being the fundamental 
Heisenberg back-action of measurement with a QPC |8]. 

Master equation approach 

As discussed in the Methods, the theory calculation involves two initial steps: 
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1. We first describe the generation of non-equilibrium "hot" acoustic phonons by the QPC 
charge fluctuations. We do this by calculating the Keldysh Green functions of the acoustic 
phonons to first order in the electron-phonon coupling to the QPC. This coupling Hamilto- 
nian takes the form: 

i?int,qpc = E ^ Q (Vf + (10) 

The "heating" correction to the phonon Green functions can be expressed in terms of the 
charge noise spectrum Sqq[(jj] given above. Here R denotes the position of the QPC with 
respect to the midpoint between the two quantum dots (see Fig. 1). Similar to our treatment 
of the dots, the spatial extent tq of the QPC charge distribution serves as a high-energy cutoff 
in the above interaction; we take this scale to be small enough that it plays no significant 
role (i.e. the QPC voltage instead provides the relevant cutoff). 

2. We next calculate Golden rule rates r^,r^ describing transitions between the DQD states 
1^) and |e) via the DQD-phonon interaction given in Eq. 1 of the main text; this is done 
using the "heated" phonon Green functions computed above. 

Finally, we incorporate the rates r^,r^ into a master equation describing the probabilities 
of the states \g), \e) (which have a total of 3 DQD electrons), as well as the two-electron state 
(1, 1) (denoted |2)) and the four-electron state (2, 2) (denoted |4)). This is similar to the approach 
outlined in Ref. [llj. This master equation describes the electrostatic blocking mechanism depicted 
in Fig. 4 of the main text. In addition to the phonon-assisted rates, the master equation also 
involves rates describing incoherent lead tunneling to and from the DQD. A slow incoherent rate 
Fgiow (involving lead tunneling from the right) describes transitions from state |2) to |^), and a fast 
incoherent rate Ffast (involving tunneling to the left lead) describes transitions from state |e) to 
These rates (which are set by the tunnel coupling to the leads) also determine the incoherent 
rates describing transitions from |e), \g) to the (2,2) state. 

We are interested in the experimentally relevant limit where Ffast ^ ^siowi corresponding to the 
conditions underlying the blocking mechanism. In this limit, and for gate voltages far from the 
charging-hnes for the (2, 2) state, the stationary probabihty Pg to be in the DQD ground state \g) 
takes the simple form: 

. (11) 
Importantly, the fast rate Ffast does not enter the limiting equations when Ffast ^ Tsiow 

In 

addition, it is straightforward to show that the relaxation rate connecting states |e) and \g) does 
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not affect this result so long as <C Ffast- The magnitude of the ground state occupation therefore 
depends only on the relative magnitude of the slow refilling rates and the coherent excitation rate. 

Finally, we can use the master equation to calculate the derivative dneff/dV/,, where Vl is the 
gate voltage used in the experiment to extract the differential transconductance (see Fig. 3 of the 
main text), and 

^eff = {nn) +e{nL) (12) 

is the effective charge sensed by the QPC. This quantity is proportional to the experimentally 
measured differential transconductance d/qpc/dV/,, thus allowing a comparison between theory 
and experiment. The parameter e — 0.4 is determined experimentally from the QPC's relative 
sensitivity to charge addition to the L dot versus the R dot. 

Numerical and experimental parameters 

The theoretical plots in the main text take the distance between QDs in Fig. [8] to be = 
If/, — r/^l ^ 235 nm; this is in reasonable agreement with estimates made from SEM images of the 
device, and also yields a spacing between back-action-excitation lines in the stability diagram that 
match experiment. Based on estimates from device images, we take the DQD-QPC separation 
to he R 500 nm. Electrostatic energies used in the theory are obtained from experimentally 
measured charging diagrams. We find charging energies of the left and right dots to be Ec^l — 2.9 
meV and Ec^r — 2.7 meV , respectively, while the interdot charging energy is Ea — 0.5 meV. We 
also use a value of interdot tunnel coupling tc ^ eV (being half of the energy splitting between 
|e) and \g) for 5 = 0) that is extracted from measured stability diagrams. Finally, as already 
discussed, the theoretical calculations take the spatial extent of dot and QPC charge distributions 
to be small enough that they play no role (rg = 2 nm). 

The DQD experiments described in this article employ a small QPC transmission T ^ 0.0028 
in order to avoid back-action due to shot noise and heating effects due to large QPC powers. 
The change in transmission associated with a change in the DQD charge state is measured to 
be AT ^ 0.00043. Via a simple calculation using a screened Coulomb potential, one can show 
that moving the excess electron from the left dot to the right results in a potential change of 
/SU ^ 25 fiY at the position of the QPC, which is consistent with experimental estimates. The 
QPC biases employed in our measurements are on the order of Vqpc ^ 1 mV. This procedure is 
used as a convenient way to determine the slope of the QPC's transmission curve as a function of 



23 



local potential changes. 

For the incoherent rates Ffast and Fgiow connecting the undriven DQD to the leads, and entering 
our master equation calculation, only rough estimates are available from our experiments. For our 
numerics, we take Fgiow = 10 kHz, and set Ffast — 1 GHz such that Ffast ^ Tsiow We stress that 
in this limit, the specific value of Ffast plays no role. The value of the slow rate Fsiow only serves to 



determine the overall magnitude of the back-action-induced probability oscillations, cf. Eq. (11). 
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